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TECHNICAL NOTE 2837 


CORRECTIONS FOR DRAG, LIFT, AND MCMENT OF AN AXIALLY 
SYMMETRICAL BODY PLACED IN A SUPERSONIC TUNNEL 
HAVING A TWO-DIMENSIONAL PRESSURE GRADIENT 
By I. J. Kolodner, F. Reiche, and H. F. Ludloff 


SUMMARY 


The corrections for drag, lift, and moment are derived for an 
axially symmetrical "body placed in the test section of a supersonic tun- 
nel, on the assumption that the test section is characterized hy a two- 
dimensional pressure field originating from construction flaws. Althou^ 
relatively simple longitudinal and transverse pressure gradients are 
assumed, the analytical treatment hecomes rather difficult because of the 
difference in symmetry between the body and the basic flow field. 

Assuming irrotational conditions, the velocity potential of the flow 
around. the body is expanded in a threefold manner: (1) In powers of the 

thickness parameter of the body €, (2) in powers of a parameter b 
characterizing the inhomogeneity of the basic flow field, and (3) as a 
Fourier series in the azimuth 9 around the body axis. Each expansion 
is taken into account not further than up to the second term. 

Upon substitution of this potential series, the nonlinear equation 
of motion and the boundary condition on the body surface are split into 
a set of linearized boundary- value problems which can be solved analyti- 
cally. The mathematical techniques used for the solution are explained 
in appendixes. 

Assuming the two-dimensional pressure field, the drag, lift, and 
moment corrections for arbitrary body shapes are obtained in closed ana- 
lytic form. The physical meeming of the results and their validity are 
discussed. 


INTRODUCTION 


Consider an axially symmetrical body placed in the test section of 
a supersonic tunnel. The test section. Instead of providing uniform 
flow, may be characterized by a two-dimensional pressure field. In 
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general, such a pressure field will consist of a longitudinal as well 
as a transverse pressure gradient, producing a stream-angle variation 
along the tunnel axis. 

The difference in symmetry between the body and the field engen- 
ders considerable difficulty in the analytic treatment of the problem, 
involving a Fourier expansion of the disturbance potential of the body. 
Therefore, a relatively sinqple, linear gradient is assumed: On the axis, 

the horizontal component of velocity may equal the original velocity of 
the uniform stream Uq, but the transverse gradient may produce a verti- 
cal velocity con^jonent 

q.y = 'bx (1) 

yielding the desired stream-angle variation. Here b is a small param- 
eter which characterizes the first-order deviation of the actual basic 
flow from the uniform field. 

This work was done at Hew York; University under the sponsorship and 
with the financial assistance of the Habional Advisory Committee for 
Aeronautics. 


SYMBOLS 

b parameter characterizing inhomogeneity of pressure field 

b = bZ/Co 

c local sound speed 

Cq sound speed corrresponding to 

D drag 

f (x) body profile function 

k(x) = 

L lift 

Z body length 


M 


moment 
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Mo = 

P 

Po 

+ 

Uo 

x,T,e 

x,y,z 


o 

local static pressure 

static pressure at velocity Uq 

v^ + 

original velocity of basic field 

axial, radial, and circumferential 

« 

cylindrical coordinates 
rectangular coordinates 
local angle of attack 


velocity con^jonents 


p = - 1 

' y adiabatic exponent 

e body thickness parameter 

p local density 

Pq density at velocity Uq 

cp potential function 


TWO-DIMENSIONAL BASIC FLOW FIELD 


The two-dimensional velocity pontential cp of the basic flow field 
satisfies the well-known potential equation: 


(c 2 - - 293 t<Py<P^ + (c 2 - <Py2)Tyy = 0 


( 2 ) 
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where c is the local velocity of sound. The Bernoulli eq.uation yields 



where 

2 2 2 

= cp/ + V 

Cq is the value of c when q = 8Jid 7 is the adiabatic 
exponent . 

Now assume that the basic potential may be represented as a foimal 
power series in b: 


cp = cp^ + bcp^^ + b^®^ + . . . 


( 4 ) 


This procedure will permit determination of the functions cp^®, cp®^, 

and so forth. It is immateidal whether this procedure will lead to a 
convergent power series, since the results may be considered as satis- 
factory even if the series has only an asymptotic character. 

Choose 


cpOO 


UqX 


and require that for y=0 (i.e., on the x-axis) 

'Px = 

<Py -ta 

(see "introduction"). This implies that, for y = 0 , 



cp^Ol = 0 


cp/1 = X 

cp^o^ = = 0 




for n > 1 


( 5 ) 


( 6 ) 


( 7 ) 
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Inserting equation (4) into equation (2) and requiring that this he an 
identity in b, it is found that, up to second-order terns in h. 


01 _ 01 = 0 
^xx ^yy 


(8) 


where 


^ - 1 = - 1 
^o 


and 


O2._»?o 


yy c 


o i- 


(7 + 1 ) + (7 - i)p‘^ 


fPx^^xx°^ + (9) 


The general solution of equation (8) is 

= f(x -pQy) + g(x + p^y) 

where f and g are arbitrary functions. 

In order that conditions ( 7 ) he satisfied, one must have 


f(x) + g«(x) = 0 
-P (f» (x) - g('x;)] = X 

This in^jlie’a that 

f(x) = -g(x) 

from which it follows that 


= xy 


( 10 ) 
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Inserting this result into equation (9) gives 


02 

P 'Pxx 


cpyy02 = -2 ^ X 


(9a) 


Mn 3 , . 

Now^ ^ — X is a particular solution of equation (9a). Hence, the 

general solution is 


cp02 = _ + F(x - py) + G(x + Py) 

3P^Cq 

» 

•where F and G are arbitrary functions. 

Proceeding as before to satisfy conditions (7)^ oae easily finds 

that 


, F(x) . G(x) - 1§— x3 


Substitution yields for the slii 5 )le form 


02 Mq 2 
cp = ^ xy'^ 
Co 


( 11 ) 


From equations (i<-), (lO), and (11), 'the potential of the basic flow 
becomes: 


9 


= UqX + bxy + ^ xy^ b^ + , . 


basic 


( 12 ) 


Hence : 


9, 


= Uq + by + ^ ySijS + . . 




M. 


2 

cp = bx + 2 — xyb + . . . 
y Cq 


(13) 
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The equation of the streamlines, 
turns out to he: 


in the present approximation. 


y - 


»oCo 



1 


^ 0^0 




+ Constant 


(14) 


From the Bernoulli equation (equation (3))^ in comhination with the 
adiabatic relations 



(15) 


where p is the density, one obtains for the excess pressure Ap = p - 
the expression: 







(16) 


Using equations (13) ^ this yields, up to second-order terms in b. 




(“o‘ 



(17a) 


and on the x-axis. 


^ = -1 x2b2 

Po 2 


(17b) 
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Therefore the pressuxe gradient of the basic flow has the two components: 


Vy.= (Ap)x = 


p„ = (Ap)„ = -Po^o^ - po(^o^ + 




(18) 


In the manner indicated above, all the hi^er-order coefficients in 
the e2q)ansion for cp could be determined successively. (They turn out. 
Indeed, to be polynomials in x and y. ) These coefficients are, how- 
ever, not needed in the present work. It is understood that expres- 
sions (12), (13), and (17) for the basic potential, the corresponding 
velocity components, and the excess pressure are valid only in a 
restricted region of the working section around the disturbing body. 


AXIALLY SYMMETEICAL SLEMDER BODY IN TWO-DIMENSIONAL 
BASIC FLOW - THE DISTURBANCE POTENTIAL 
Differential Equations 


It is assumed that the total velocity potential cp can be approxi- 
mately written as 

<p = ( 9 °° + bpOl + + (,V° + ■ (19) 

This assunption is by no means trivial nor arbitrary. It is justified 
by the fact that in this way it is possible to satisfy the bounciary con- 
ditions on the body and the "characteristic condition" (see appendix A) 
in such a way that the error committed- in cp' as given by equation (I9) 
is small ctmpared with the teims written down. The method adopted here 
is sucdi that it can be extended (stepwise) to better and better approxi- 
mations. It is not obvious, however, what the dependence of the next 
term on b and e is. The method, if pushed further, leads to a 
formal series in b, €, and logg g, with, at best, an asymptotic char- 
acter. For further elucidation, see the section "Boundary Conditions." 

The terms in the first parentheses in equation (19) represent the 
basic flow potential as discussed in the section "Two-Dimensional Basic 
Flow Field," while those in the second parentheses represent the disturb- 
ance potential produced by the body, and g is the thicikness parameter 
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of the body whose equation, written in cylindrical coordinates (r, 0, 

x), is assumed to he 

r = ef (x) 0 ^ X ^ Z 

f(0) = f(Z) = 0 


The first term in the disturbance potential is the well-known term of 
the linearized theoiy. The second term represents the influence of the 
pressure gradient (of the basic flow) on the disturbance potential. In 
the present work, higher- older approximations will not be considered. 

An important and rather restrictive assumption has been tacitly 
made by using expansion (19) for cp, namely, that 

logg e « b « 6 (20) 


(In order to conqpare different orders of magnitude, the dimension- 
less quantity b = bZ/CQ is introduced. In the equations in the text, 

the dimensional parameter b is kept, and inequality (20) is taken into 
consideration in a suitable manner. ) 

That expression (20) follows from equation (I9) can be seen from 
the following argumentation; If one sets 6=0, the potential cp must 
reduce to the basic potential, since a needlelike object produces no 
disturbance ln_a three-dimensional flow. On the other hand, if one 
sets b = 0 (l.e., assumes that the basic flow is uniform), cp must 

reduce to the disturbance potential of a body moving in uniform flow. 

It is well-known (see reference l) that the second term of this poten- 

tial is of the order loge 6. This term has been omitted in equa- 
tion (19) as small congjared with a term in e^b, and this is permissible 
only if b » 6^ l°ge other hand, the term in b^ in the basic 

potential has been omitted as small compared with all terms written down, 

in particular, TOth the smallest term be^. This implies that b « e. 
Without this restrictive assimrption (20) the computational work would 
become prohibitive. 

The cylindrical coordinates x, r, and 9 are introduced in such 
a way that the plane of the -basic flow is the plane 0 = 0. Then, 
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putting y = r cos 0 and z = r sin 0, the potential equation for 
q)(x, r,0) becomes 

f99 ' ^ 'PrWrS ' 

^ fx'Pel’xe - ^xVxr ■" P * ^ ‘Pe^)r "fr ' ° 

Q 0*0^1 

Here, c is again given by equation (3), where now q = ^ 

Inserting expansion (I 9 ) into potential equation (21), and ordering terns 
in powers of b and 6, one obtains in orders b and b^ equations (8-) 
and (9)j in order 6^, 


„2 10 10 1 ^ 10 „ 

Bcp -CD -=-Cp =0 

^xx ^rr r ^r 


( 22 ) 


in order e^b. 


2 11 11 1 11 1 11 
-) _ (D _ m ^ m 

XX 


P cp_ - - F 'Pr - ± 'P 00 


M. 


[ (7 + 1 ) + (7 - l)p^ 

rcp^lO + 2xtp^l0 + 2cp^l0 V 

u- 

V. 


cos 0 ( 23 ) 


Equation (22) is the well-known equation for the linearized poten- 
tial cp^®. The term 1ms been omitted since, as is well-known, 

J.2 00 

the linearized potential does not depend on 0. Equation (23) for the 

"interaction potential" cp^^ is Inhomogeneous and, since its ri^t-hand 

side involves 0 explicitly,* cp^^ certainly depends on 0. Expanding 

cp^^ into a Fourier series (since the problem is symmetric in 0, cp^^ 
is obviously a cosine Fourier series) there ip obtained: 
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= y cp^^^(x,r)cos n© (24) 

n =0 


Substituting equation (24) in equation (23), it is found that the Fourier 
coefficients satisfy the differential equations 


^ ° 1 ( 25 ) 


»% - <p - I <p + i ,p 11 s(x,r) (26) 

-i-xx -‘-IT ^ -*-r r^ ^o 


where 


S( 


x,r) = j7y + 1 ) + (7 - l)3^rcpjQjlO + 2 q >^10 + 2 xrp 2 ^ 1 ° 


(27) 


It will he shown in the next section that the boundary conditions can be 
satisfied only by making - 0 for n / 1. Hence the interaction 

potential reduces to 


cpll = cp^ll(x,r)cos 9 


( 28 ) 


Boundary Conditions 

The flow velocity must be tangential to the body at its boundary, 

hence 

= ef’(x)cpjj (29) 

for 


r = ef(x) 

Inserting expansion (19) in equation (29) and taking account of eqiia- 
tions (5), (lO)’, and. ( 11 ), the boiindary condition becomes: 



12 


MCA' TN 2837 


p / M_ p \ O ICS 2 — 11 

bx cos 9 + D 12 — xr cos ») + ecp-f. + be ^ cpu cos n0 + 
^ ' n=0 

/ 2 2 ? 2 10 
ef‘(x)jUQ + br cos 0 + b ^ r cos 0 + e q>^ + 



n=0 



cos n0 + 


(29a) 


for 


r = ef (x) 


To be able to order the terms of equation (29a) in an appropriate 
manner, it has to be known how depend 

on 6 for r = ef (x) . 

This knowledge is gained by considering solutions of the differ- 
ential equation 

P^xx - ^rr - ? ^r + 4 ^ =0 (30) 

r 

which vanish to the left of the cone r = x/3 and tend to zero as one 
approaches Infinity along a characteristic r = x/p + Constant. It is 
shown in appendix B that such solutions are given by 


2^r“n' UQ 


- I) + - 1)^ - -I- jjx - O - V(x - if- - 


« (31) 


\/(x - if - 
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where g(| ) is an arbitrary, n times differentiable function which is 
identically zero for ^ < 0, and that, as r — > 0, 



It may now be shown how the boundary conditions for cp^^, and. 

so forth can be determined successively in a unique way. First, 

satisfies equation (30) with n = 0; hence cpj,^*^ = ^ ^ 

and is of the order e for r = ef(x). Matching terms in e 

in equation (29a), one now gets the condition 

lim rcpj.^® =Uof(x)f»(x) (33) 

r — ^ Q 

This condition deterniines cp^® uniquely (see the section "Solution of 
Boundary-Value Problems"). Next consider which also satisfies 

equation (30) with n = 0- Using equations (32), the term be^cp_ 0(bc) 

Up 

as r = ef(x)j since there are no other terms of this order in equa- 
tion (29a), -one must set = 0, which of course is a solution of 

equation (25). Now, the terms for n > 1 are of the order 

be^”^ and become infinite as e — ^ 0, unless V = 0 for all values 

of n greater than 1. This proves the statement made at the end of the 

section "Differential Equations." The potential cp^^ does not satisfy 

equation (30) but is a solution of an equation with the same left-hand side 

with a nonvanishing right-hand side. Thus cp^^^ is obtained by adding 

to an expression of the type of equation (31) a particular solution of 
equation (26). It is shown in the section "Solution of Boundary-Value 
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Problems" that it is possible to find such a particular solution with 
the r derivative of hi^er order than ^ for r — ^ 0. Hence the 

singularity of cp^ is that of ijr for n = 1, that is, cp, 

ir ir 

as r — ^ 0, and the term 

be^cp, cos 0 = ^ — 5- lim cos 0 + 

Terms vanishing to higher order 



The only term of this order in equation (29a) is bx cos 0 and so one 
gets the condition: 


lim 

r — ^ 0 


= -xf^(x) 


(34) 


Grouping the matched and unmatched terms in equation (29a) and 
substituting r = ef(x) wherever r appears explicitly, equation (29a) 
now reads; 


1 

f(x) 


xi^(x) 


lim (rcpr^^) - Uof(x)f’(x) 


e + 


f^(x) 


lim 

r — > 0 




b cos 0 -f — xf(x)b^6 + — xf(x)b^e cos 20 - f’(x)f(x)be^ cos 0 - 




2c 


2- f* (x)f^(x)b^e^ - ^ f’(x)f^(x)b^e^ cos 20 + 


2 c 


10 
e cp 
^r 



+ 





11 


e^cp^lOf > (x) - e3bf « (x)cpj^^ 


COB 0 + 


0 


(35) 
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The first two brackets in equation (35) are zero^ by virtue of the 
iii 5 )osed boundary conditions (33) and (3^). All the unmatched teims are 
of hi^er order in g and b and can be disregarded. In particular,, 
the terms in the last two sets of brackets are of order hi^er than g 
and b, respectively, whereas the last two terms are of the order 

log 6 and 6^, respectively, as may be seen by equations ( 32 ). 

These terms could be matched only by considering higher approximations 
to cp. 

The preceding discussion makes clear how to proceed with hi^ex- 
order approximations . The various terms in the expansion all satisfy 
equations of the type of equation ( 30 ) with nonhomogeneous terms. At 
each step only a finite number of Fourier coefficients will be con- 
sidered. The singular behavior of these Fourier coefficients becomes 
worse with the order of the coefficients, but this is in turn compen- 
sated by an appropriately large exponent of 6 in the expansion, thus 
making the coefficients appropriately small at the boundary of the body. 


Solution of Boundary-Value Problems 

The solution of equation (22) with boundary condition (33) is the 
well-known expression for the linearized disturbance potential of the 
Karman-Moore problem, -vdiich is used herein in the fom given by Courant 
and Friedrichs (see reference 2) 


where 




k'(g) dl 
\/(x - I) - pV 


k(x) = i f^(x) 


(36) 


This is exactly the solution ijr of equation (31) with n = 0 and 
g(x) = i Uok’(x). Clearly, using equation (B6b) of appendix B, 

lim = 2g(x) = UQf(x)f ' (x) . It is assumed that: 


f(0) = f(Z) = 0 


(37) 
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Tills Implies that 


k(0) = k(2) = 0 
k*(0) = k*(l) = 



(37a) 


Using foimulas (Bl) and (b 4) of appendix B^ upon setting n = 0 and 
g(x) = |-Uok’(x), it follows that: 


<P “ = 




r k"(i) ae 
° iix - if - pV 


k" (k - pry) dT 

\/t2 - 1 








10 = ^ 


^ (x - l)k(l) d| 

= II B Tk(x - Prr) dT 

° \/t2 - 1 


(38) 


To solve equation (26) for subject to boundary condition (34), 

one must first find a particular solution of equation (26) and add to it 
an appropriate solution of equation ( 26 ) with the rl^t-hand side equal 
to zero (homogeneous case), such that the sum satisfies th,e required 
boundary condition, equation (34). It is verified in appendix C that a 
particular solution of equation (26) is given by 



_Mo 


rcplO + xTcp^l® + 


Mo (7 ^30 


4P' 


(39) 
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It is observed that: 


lim r‘ 
r — > 0 




m 


This follows, since cplO = o(logg r ), qp^lO = O(logg r), cp^lO = 

= ^(r)^ virtue of equations (32) with n = 0. 

/ -1 1\** 

Hence the appropriate solution of the homogeneous equation I y nrust 

satisfy the same boundary condition as that imposed on itself; 

/ 

(qpj^ 1 is now immediately fotmd, since it satisfies equation (30) with 
n= 1, and 

lim = -xf^(x) = - 2xk(x) 

By equation (B6a), appendix B, the desired solution is obtained from 
formula (31) with n = 1 and g(x) = -2xk(x). Hence: 


L ll\** =2 (x - O |k(|) + |k«(g J dj 

\l J ^ T Jq r 




(41) 


Using equations (39)^ (4l), (36), and (38), there is finally obtained 


V" ■ - (V")“ 


^2 [k'(D + Xk*'(5)] dl 


■f 


^(x - D" - 


p2 2 

) - P r 

^ „x-pr (x _ 0[^(0 + 2gk'(5) +Pr2k''(|[| dS 


iL 




(x - - pV 


= -Mo ^ 


2 pVP^ [k’(x - Prr) + xk"(x - PrrjJ dr 




x/Pr t[^(x - Prr) + 2(x - Prr)k' (x - prf) + Pr^k”(x - Prx^ dx 




i 


( 42 ) 


x^= - 1 
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vhere 

Mo^(7 + 1) 

P = — 


Differentiating, one gets: 




2 I— /P [i''(x - prr) + 3k***(x - Prr^ dr 




^x/P Tf4k'(z - Prr) + 2(x - PrT)k"(x - Prr) + Pr^k” ' (x - Prr) dr 
P / _i= =: + 

\/t 2 - 1 

(P - «o^) 


XT V* 


1^x2 . pV 


k”(0) 


(43) 


px/pr -^[^'(x - Prr) + xk''(x - prT)| - prr[k**(x - pir) + xk’"(x - PitTIV dr 

L ^-7= ^ ^ 

px/Pr -|2PrTk"(x - Prr) - Pr2[kk'(x - prr) + 2(x - prr)k”(x - prr) + Pr^’"(x - Prr)]j^ dr 


(2 - «o^) 


\/x2 - p 2i2 


k"(0) 


(44) 


and obviously: 


cos 9 


qp^ll = cos 9 

cpgll = sin 9 


> 


(^ 5 ) 
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COMPUTATION OF VELOCITIES, PEESSUTRES, AND FORCES 
Velocity Con5)oneiit3 on Surface of Body 


Denote by u, v, and w the axial, radial, and circumferential 
velocity components, respectively, and split each one into three parts, 
due to the basic potential, the first-order disturbance potential, and 
the interaction potential. The first of these are computed immediately 
from equation (12), whereas the computations involved in evaluating the 
others are carried out in appendix D. 

In the following relations it is important to have all terms arranged 
in the proper oixier of magnitude. In view of relation (20) 


S logg 


e « b « e 


one has 

1 » € » S » logg € » logg € » log^ 

6b »b^ »b6^ logg 6 » 6^ logg^e » e** log^ € 

6^ logg 6 » 6^ » 


_ 2 

6 » be 


» b^€ » 


> ( 46 ) 


J 


Up- to teirms in be^, one now has, by differentiating equation (12), 


'^basic = Uq + b6)/2k(x) cos 0 + . . . 


^basic = "bx cos 0 + .. . . 


y (47) 


^basic “ sin 0 + . . 
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Using formulas (D8) and (D9) of appendix D: 


loge - .%o|r 


fiSr' 


+ Terms of order logg e 


=2(p 10 ^ + g3 eUoP^l^p^ k* • ’ (x) - 

e V 2 


g3u |y "" S»««(x - g) - k«««(xj dn ^ k««(0) ^ 


> m 


k*“(x) 






€2(p^l0 = 0 


Using formulas (DIO) to(Dl2) of appendix D: 


1)6^, 11 cos 0 = "be 

-Lv 


^IS(7) + 

V5w J 


COS 6 +., . 


11 cos 6 = -bx cos d + be^ l°Se ^ [(3Mo^ - 2)k'(x) + 

(2Mq2 - l)3i:”(x^ cos e + bc2Q(x) cos 9 + . 

= -bx sin 0 - b6^ l°Se ^[(3Mo^ - 2)k*(x) + 

(aio^ - l)xk"(x^ sin 0 - be^*(x) sin 9 + .. . 


• y (49) 
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Frcm equations (47), (48), and (49) one obtains: 
~ '^basic '^disturbance 


= Uq + 6^ logg eUok“(x) - e%o 




+ 2 be 

\/aw + f W 

logg e 

]/2k(xl 


k''(x - g) - k’»(x) 


cos 6 + 


dcf + 


Terms of order loge ^ 

At present all one needs to know about w and v is that ' 

w = - 2 bx sin 0 + ^ ® "*■ T®^™s of higher order 


( 50 ) 


(51) 


and that 


V = e + “1^ ^ ^°Se ^ + Terms of hi^er order 


V^x) 


It is i mma terial idiether one knows what the above braces actually 
represent. 

Now consider again boundary condition (29a) 

cpj, = ef ’ (x)q) 3 j for r = ef (x) 


(52) 


or 


V = ef'(x)u = e ^ u 

V 2 k(x) 


Substituting u and v obtained in equations ( 50 ) And (52), it is now 
observed that the unma tched term of lowest order is proportional to 

o 

^ e. It may be easily verified by use of the method qf the sec- 

tion "Boundary Conditions" that introducing a new term in the potential. 
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to cancel this teim in the boundary condition, will result in the addi- 

I4- p 

tion to u of a term of order e logg‘^ 6 ; hi^er-order corrections will 
introduce in u terms of still hi^er order- Consequently, formula ( 50 ) 
for u is certainly correct to the order loge ^ (see rela- 

tions (46)). The sa 3 ne' ca n not he said about formula (52) for v and 
formula (51) for w. In these equations, hi^er-order potentials will 

introduce corrections of order loge so that the method yields 

expressions for v and w correct, on the body, only to order logg e. 
This is already sufficient to carry out the computation of the pressure 

field near the axis, correct to the order be logg e (see next section). 
On the other hand, one can easily find a better formula for v, on the 
body, than the previous one. Indeed, using the boundary condition and the 
present knowledge about u, the expression 


V 


6 

\/ 2 k(x) 


u 


(53) 


is correct to the order be^ logg e, hence definitely to the order 


be' 


The quantity 


S 




» 


where 


q^ = + 


v2 + 


w2 


(54) 


is of importance in the evaluation of pressure on the body, drag, and 
lift. Since the expansion for u begins with a term of order unity, 
this quantity may be conputed at best to an order not hi^er than or 

equal to logg^g. ' Using equations (50), (53), and ( 5 I) one easily 
finds: 
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S = logg 6Uo%”(x) + e%, 


•■{/: 


k«^'(x - 0) - k»«(x) ^ 


I — [k‘(xyi^ 

k’'(x) log^ ^J^-r - -^= =^> - cos 0 

Plfk(x) 4k(x) 


\/2k(xJ 


xk’ (x) 


\/2k(x) 


- b^(2x^ sin^0) + 


(55) 


correct to order be^ log^. e. 


Pressure, Drag, Lift, and Pitching Moment 
Using formula (I6) for the excess pressure, one finds 


AP = PqS 


1 


P = Po + PoS 


(56) 


correct to the order be^ lege 

The drag is the x cem^jonent of the pressure force integrated over 

the body surface. The surface element of the body is where 0 

cos a 

is the azimuth angle and a, the local angle of attack. Then the element 
of the drag force is 


dD = pr tan a dx d0 


Since 


pr tan a = pe^f(x)f'(x) = e^k'(x) 


one gets: 


D 



pk’(x) d0 


(57) 
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As 



dx = 0 



k'(x)k’*(x) dx = 0 



cos 6 dS = 0 


only terms in and in of S, equation (55)^ contribute to tbe 

dreig. Thus, one obtains 


D = D]_ + Dg 


(58) 


Here 


Di = 


-2npoUo26^ 




k' *(x)k* ’(cr) logg (x - a) da 


(59) 


is tbe well-known formula for tbe drag in tbe linearized theory (see 
reference 2, section 153; for a simpler computation, see appendix E) and 

n Z p2it 

®2 “ ^ / 2x^k’(x) sin^6 dJ9 



xk(x) dx 


(60) 


where ©2 represents the effect of the inhomogeneity of the flow. 
Observing that 
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is the volume element of the body, one may vrite: 

D 2 = 2Po5Vh^ 


( 61 ) 


where x is the distance of the nose of the body from its center of 

I 

.2 / ^ volume of the body. Formula (58) 


mass and V = 2it6‘ 


' f k(x) 


for the drag is correct to the order bc^ l°Se 

The lift is the y component of the pressure forces integrated 
over the body surface. The element of lift is 

dL = -pr cos 0 dx d0 

= -ep\/2k(x) cos 6 dx d0 

All terms of p contribute nothing to the lift, except for the term 
in 6b (see equations (55) and ( 56 )). Thus 

cos^ d0 

= PoUoVb ( 62 ) 


an expression which is correct to the order be'^ loge 

The element of the pitching moment about the nose is: 

dM = p cos 0 rx dx d0 

= 6p^21c(x)x cos 0 dx d0 


Again, all terms of p contribute nothing to the moment, except for the 
term eb (see eqiiations (55) and ( 56 )). Hence 


M = -kpgUgbe^ 





cos^0 d0 = 0 ( 63 ) 


an expression which is correct to the order bs'^ l°ge 
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Validity of Results - Physical Interpretation 

As was stated at various points in the text, the derivations do not 
yield results which are all correct to the same order of magnitude. For 
clarity, the main conclusions are listed again. 

The main assimption on orders of magnitude is that 

logg e « h « e 

The potential and velocity coii5)onents away from the hody are determined 

p 

to the order he . They contain terms of orders 

2 ? 2 
1 » h » e » b » be 

On the body they are determined as follows; 

(1) The u coa^onent of velocity is determined to the order 

p 

be logg e and contains terms of orders 

1 » logg 6 » » be 


( see formula ( 50 ) ) • 

(2) The V con^jonent of velocity is determined to the order 
be-' logg s and contains terms of orders 

e » e logg e » 6^ » be^ 

( see formula . ( 53 ) ) • 

(3) The w con^jonent of velocity is determined to the order b and 
contains only a term in b (see formiila (5l)). 

(U) The pressure is determined to the order be^ loge ^ contains 

terms of orders 

1 » e^ logg e » e^ » be » b^ 

(see formulas (56) and (55))* 
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The drag Is determined to the order he^ loge e and contains terms 
of orders » h^e^(see formula ( 58 ))- Both the lift and the pitching 
moment are determined to the order hg^ logg e, the former containing 

p 

the only term of order he and the latter being zero to this order. 

Of the two drag terns, the one proportiohal to is the wave 

drag of the body, in the representation by Courant and Friedrichs. 

One mi^t anticipate from general physical considerations that a 
"horizontal-buoyancy" term will occur which should be equal to the 
product of volume and pressure gradient in the x— direction. From equa- 
tions (18) it can be seen that for the assumed basic field, the x gra- 
dient of the pressure, up to order b, is zero. Therefore no drag term 

p 

of the order be may be expected. Note that in order to con^jute the 
drag, the expression for the pressure has to be multiplied by since 

the surface element is proportional e^. 

The second drag term occurring in equation (58), of order b^e^, 
can be seen to originate from two different terms in the velocity poten- 
tial. One contribution comes from the so-called interaction poten- 
tial qp^, having as a factor be^; Itself is shown to be propor- 

tional to l/r. The corresponding pressure term turns out to be propor- 
tional to (qPr^^) and is thus proportional to b^ near the axis, so 

p 2 

that the drag term becomes of order b e . The other contribution origi- 
nates from the potential coefficients qfP^ and qfP^ of the basic field, 

which are proportional to b and b^. Since the basic field has no 
singularity on the axis, the corresponding velocities and pressures do 
not change their order, as the body surface is approached. Therefore, 

also, these- terms will make contributions to the drag of order b^e^. 

A drag term of' order b's 3 -which could be anticipated from contri- 
butions of qfP^ (proportional to b) and of qp^ (proportional to be^) 
can be shown to vanish when the pressure distribution is integrated over 
the body surface. 

As seen from eqtiation (62), the lift consists only of one term which 
is proportional to bc^. The origin of this term can be traced back to 

contributions from qfP^ as well as qp^^j q[P^ leads to a pressure term 
which is proportional to y so that the whole expression for the pres- 
sure on the surface becomes proportional to be and the lift expression 
becomes proportional to be^^ qp^^ leadn to a pressure teim which is 
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proportional to thus altogether of order he, hence again 

p 

a lift term of order he results. 

One might expect that, lift terms of order he, he^, and e^ occur, 
in analogy to the corresponding result for two-dimensional bodies. But 
it can he shown that the respective expressions vanish when the corre- 
sponding pressure terms are integrated over the hody surface. 

The moment about the nose of the hody has been shown to vanish, up 
to order he 3 logg ,e. This peculiar result, which holds for an arbitrary 
hody contour, must he attributed to the very special basic field that has 
been chosen. 


New York University 

New York, N. Y., November 5, 1951 
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APPENDIX A 

CHAEACTERISTIC CONDITION 


It must, "be postulated that the total disturhance potential 

+ he^cp^ (A1) 


vanishes at the "characteristic surface," the equation of which may he 
written in the form 

r = fl(x) = + hD^^^(x) + . . . (A2) 

where 

= x/Po (A3) 

and Pq^ = - 1. Inserting equation (A2) into disturhance poten- 

tial (A1) and using a Taylor expansion, there is obtained, neglecting 
higher- order terms. 



Hence; 


and 


(plO 


xjfl^*^^(x^ = 0 


(A4) 


(x^ + qp^ x;D^*^^(x^ = 0 (A5) 


Using the notation of appendix D, 

cpl° = -Uolo(k') 

= UoPoIi(k") 


(a6) 


qpr^ = cos 
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or 


= cos e|^Mo^rIo(k*) - xiMo%o(k") + 2PoIl(k) + 

2PoXll(k') + PPor^Il(k* ') - 2Po^rl2(k'^ (A?) 

( S66 equation (^2) .) It is easy to prove that all the integrals In(s) 
approach zero when x approaches Pq^* ^ present case^ only the 
special values n = 0, 1, and 2 are needed. Now 


io(s) = f 


Il(g) = 


12(g) = 


Put x/Por =1+6. Then 


/^O^ g(x - PqTT) 

1 - 1 

/Po^ Tg(x — Po^t) 


dT 


1 - 1 

. T^g(x - PqTt) 

Vt2 - 1 


dr 


> 


dr 


(A8) 


iio(rfi - r 


g(x - PoTt) 
- 1 


dr 


<G 


pl+6 

o 1 


dT 


\lr^ - 


= loge 


1 + 6 + \/(l + 6)^-1 


where G = |g(x - PoTt)! is always finite. The above expression goes 

I i max 


to zero if 6 >0. Similarly, 
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approaches zero if 


and 




T dr 
- 1 



1 


where = jxg(x - finite. Also this integral approaches 

zero if 5 — ^ 0. Hence^ all the integrals In(s) occurring in equa- 
tions (a 4) and (A5) vanish for r = x/Pq = meaning that the 

characteristic condition is fulfilled. 



32 


NACA TN 2837 


•APPENDIX B 

PPOPERTIES OF THE FUNCTION DEFINED BY EQUATION (31) 

It is shown first that equation (31) satisfies equation (30). 
differentiation of ^ is carried out easily after introducing t 
as a new integration variahle. One obtains then 




n px/pr (t + \It^ - 1) + (t - - l) 


_ r 

2%I 




g^“'(x - prr) dr 


and it follows that 


P 1 r\^ 

P -ilr - t i|r = 

^ ^xx r p2 


^n+2 


2°nJ 


px/p>r + G - - 1) 

^1 yjr'^ - 1 


-(t 2 - 


^g(n+l)(j5- _prr) 


dr 


(x + \/x^ - P^r^) + (x - \/x^ - P^r^) _ ^2^2 Jn+l), 


pn+2j^+2 


- P^r" g^“"^^(0) - 


(x + - P^^) 


- (x - Vx^ - p^r^) 




The 
X - I 


(Bl) 


- prr ) + 


n 


pn+ 2 ^n +2 


(B 2 ) 
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Now^ one may verify that the integrand on the preceding p6ige^ 



so that the right-hand side is identically zero and equation (30) is 
satisfied by ilr. Differentiating eqxiation (B1) with respect to x 
and r and transforming hack to the original variable of integration 
one obtains 






gW(o) (b4) 
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As r — ^0^ "the integraaids in equations (B3) and (bU) remain finite 
and continuous over the full range 0<|<x - Pr, except in the case 
of for n = 0. Hence 


lim P (x - d^ - ^ g(“)(0) 

> 0 ■ 


=-^ r (x - 1) 

nl Jn 


n-2 (n) 


g^^U) d^ 


= -^ g'(x) n > 0 


(B5) 


r — > 0 


lim ^ f (x - (^) dl + ^ g(’^^(0) 

d Q n. 

= , \., r u - di 

(n - 1) I Jq 

= g(x) n > 0 


(B6a) 


pX 

lim (ri^- ) = 2 / g*(5) d| + 2g(0) = 2g(x) n = 0 (B6h) 

r — > 0 ^0 


provided that g(0) = g'(0) 


. . g(‘‘-i)(o) = 0. 


It follows then that 


^X = + . • . 

^r = + . . . 


n >0 > 


(b7) 


where the dots indicate terms of hi^er order in r as r — ^ 0 . It is 
shown in appendix D that, for n = 0, order of log^ r. 

That ijr, and il^p vanish for x < pr is obvious, since the 

integrands then become identically zero. It may also be verified that 
both ■'l'’x and are proportional to as one lets r tend to 

infinity along a characteristic x = pr + Constant. 
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APPENDIX C 

VERIFICATION OF EXPANSION ( 39 ) AS A PARTICULAR 
SOLUTION OF EQUATION (26) 


It is to be verified that expansion ( 39 ) 




11 \* ^ 


10 10 ^o^^^ 2 10 

+ xr® ^ + — 

43 ^ 


(Cl) 


is a solution of the nonhomogenous equation (26) 


Uf) = - <Prr - F ’’r + ^ f 




(7 + 1 ) + (7 - 1 )P^ 


10 


+ 29^10 + 2xcp^^ 


10 


(C 2 ) 


One has 

L(r2(p^l0) = P^r^cp^^lO _ r^qp^^^lO - 5rcp^^l0 _ 3(p^l0 

= ^(P^<Pxx^° - 'Prr^° - ? " 'Prr^® " 

i - ^^P^'Pxx^® 

= -4p^r9xx^° (C3) 
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+ xrcp^^° + | = 

- ^'Prr^° - + P^^^Pxxx^° + 232r<p^^° - 

xr'Pxrr^® - 3xq^r^° " = 

F "Pr^°l - - ^‘Pxr^^ = 

-2cpj.^° - 2x<I^j.^‘^ 

since - <Prr^° “ i ° (22)). Since 

m/ = p2 + 1, 


1 n 10 ^o ( ^ '*’ ^) ^2 2.0 

:(plO + xr<p^lO + r^cp/ = 


rcplO 


xrcpj,- 


10 




(7 + 1) + P (7 - 1 ) _2^ 10 

+ — P ^ 'Pr 

4p2 


From eguaiions (C3) end (C4) it follows "ttien ■that 

^('Pl^^) + P^(7 - + 2cp^^° + 2xcp^^^ 

completing the verification. 

Expression (Cl) for was found hy guessing that equa- 

tion (C2) has a solution of the form 

< 5 > = Arcp^*^ + Br^cpj,^® + Cxrcf^^*^ 
and determining the constants A, B, and C. 
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APPENDIX D 

COMPUTATION OF VELOCITY COMPONENTS ON SUEFACE OF BODY 


In order to evaluate the velocity components at the hody, one must 
find the dependence on e. of integrals of the type 


i^n(g) = J 


T^g(x - 3rr) 


- 1 


dr = 


1 g”g(x - g) 
'pr (pr)“ - p2^2 


d0 (D1) 


for n = 0, 1, and 2, 'as r = 6f(x) and 6 — > 0. Using the identities 


= -rl^ - - 3 


(Pr)2C 


grr) - 2x(x - Prr) + 


it is seen immediately that 

" 2xlQ(xg) + Iq(x^^ 




(D2) 


In order to get that part .which does not vanish as e 
Sion for Iq is needed up to the terms in r^. 


0, the expres- 


Write E = Pr and integrate 


io(s) = r 

d-R 


g(x - g) dg 




(D3) 
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three times hy parts. Since the successive integrals of ■ ? are 



one easily finds that 



in the last part the original variahle of integration t has heen used. 
Expanding the brackets in the above integrand in inverse powers of t 
(l < T < 00 as R — > O) , one gets 

(l^ + i ) loa: 2 Rt - + -v + Terms of order or higher 

\ 2 4/ ® 4 4 to; 

The contribution of the unwritten temis of the integral is of order R^ 

and may be neglected. Keeping tenns up to R^, one now gets, after 
simplifying the integral by successive integrations by parts, 

“1 2 

io(g) ' sM * ^ e-'U) loge 2 - |[e(t» ' *s’<°>](l) ‘ 

g'WH + - “) + X ^ 
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Wow 


g(x - a) p^ g(x - o) - g(x) ^ ^ X 

/ 5 dcT = / 5 da + g(x) log^ ^ 

/r Jo 


.R 


g(x - a) - g(x) 


da 


The first integral is finite; the second integral yields 


g'(x)R - ^ g"(x) + 


Similarly, 


r£ g' *(x - g) 

4 Jt, ct 


da = 2! P S».J.(x_:_a) _-_g;_'(x] ^ 

4 Jn cr 


0 


i g"(x)R^ loge ~ + . . . 


R 


Hence, finally. 


io(g) da + 


da + 


g(x) + ^ g’ '(x) 


r g"(x)R^ 
4 


lo«ef 


-||I(0) -xg’(o3]g 


(d4) 


Replacing R = Pr = pef(x) = ep\/2k(x) , one now gets, using equations (d 4) 
and (D 2 ), 


'^oie) = -g(x) logg e + r ^ da + 


•g(x) log. 


(p®) 


+ Terms of order e^-logg e + . . . 


(D5). 
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I]^(g) = ^ - ■ P g{a) dxT + (e loge 6)pippl g'(x) - 

e6V2k(x) Jo '' 2 




(x - g) - g'(x) ^ g(o) ^ 


^ i) 


+ . • • 


(d6) 


l2^S^ = — - ct) dg - i g(x) logg 
2^-ke‘^dQ 


6 + 


1 

2 


7 


g(x - g) - g(x) 




dg + 


+ • • • 


(D7) 


Using formulas (D5 ) j (d6 ) , and (D7 ) , the integrals occurring in the 
velocity components, equations ( 38 ) eind (42) to (4-5), are evaluated hy 

mere substitution. There are obtained, as r = e J2k(x) , 


cp/0 = -Uolo|k''(x)] 

= Uok’ '(x) logg e - Uq 


y 

^ 0 


k' '(x - g) - k' '(x) 


dg + 


k' '(x) logg ^ — + Terms of order logg e + . . 


PVk(x) 


(D8) 
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= UoPIijk' '(x^] = + 6 logg 6 k" '(x) - 


eUoP' 


2 k(x) 


^2k(x) 

k* ' '(x - q) - k" '(x) ^ k' '(o) 


k' ' '(x) 




+ 


(D9) 


'Pi. 


11 _ 


-Mo^\l2k(x) e.^Io|k''(x^ + xio |k' ' ' (x)^j^ + 
43I]^[k'(x)] + 23Ij^[^"(x^ + 2PPk(x)e%l[k" ’(xj] 


/ P\ xe\j2k(x) 

(p - Mn^ V k«'(0) 

Vx^ - p^e^2k(x) 


]/2k(x) 


+ 6 log e 


(DIO) 


<Pl/^ = -Mo%5'(x 3 - Mo2xiQ|k-'(x)] + 

Mo2p^2k(x) 6l^[k"(x)] + 

Mg2p^2k(x) xel]^[k"'(x^ + 2P3\f2k(x) eli[k"(x^ - 
43^12 [k'(x[| - 23^12 [xk"(x^ - 2P3^k(x)e%2[^" '(x^ - 

(P - Mg2)xk' '( 0 ) + . . . 




- 2)k'(x) + (2Mc 


1) 


xk' ' (x) 


loge e + Q(x) + . . 

(Dll) 
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where 



Dots indicate terms which either vanish as e — > 0 or tend to zero 
faster than the last term considered. In the above five formulas the 
fact that k(0) = k'(o) = 0 was ■used (see equation (3Ta-))» 
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APPENDIX E 

DERIVATION OF FORMULA FOR DRAG IN LINEARIZED THEORY 


Formula ( 59 ) for Dj^ is the well-knbwn expression for the drag of 
the linearized theoiy. For the sake of completeness, a sin^jle deriva- 
tion of this formula is given. From equations ( 57 ) ^ ( 56 ), and ( 55 ) ^ 


Dn 






'X 1,, 


k' ' (x - 0) - k' ' (x) 


d 0 + 


k' '(x) log. 


[k’(x)]‘ 


® 3 Vk(x) , 4 k(x) 


dx 


Since 


.(x) lose ^ 5^X) 




and k'( 0 ) = k'(Z) = 0 , the above integral reduces to 


= r k'(x) r 

^0 ^0 


Now, the bracket above 
k' ' (x - 0) - k' ' (x) 


k"(x - 0) - k"(x) 


d0 + k' '(x) log„ X 


dx 


I 

^0 


d0 + k' '(x) logg X = 


+ k” (x) logg X + 


|^''(x - 0) - k''(x^ loge ^ 

px pX 

/ k'''(x - 0) logg 0 = k''(0) logg X + / k’”(x 

'->0 ^0 


- 0) log 0 d0 = 


A . r k " 
^Jo 


(x - 0) log 0 d.0 
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so, after integrating "by parts, 


I = 




loge 



k' *(x)k' '(0) logg (x - 


0 d0 


0) d0 


which is the desired result. 
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